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ABSTRACT 

We estimate the evolution of the contribution of galaxies to the cosmic background 
flux at 912A by means of a semi-analytic model of galaxy formation and evolution. 
Such a modelling has been quite successful in reproducing the optical properties of 
galaxies. We assume hereafter the high-redshift damped Lyman-a (DLA) systems to 
be the progenitors of present day galaxies, and we design a series of models which 
are consistent with the evolution of cosmic comoving emissivities in the available near 
infrared (NIR), optical, ultraviolet (UV), and far infrared (FIR) bands along with the 
evolution of the neutral hydrogen content and average metallicity of damped Lyman- 
a systems (DLA). We use these models to compute the galactic contribution to the 
Lyman-limit emissivity and background flux for ~ z < 4. We take into account the 
absorption of Lyman-limit photons by HI and dust in the interstellar medium (ISM) 
of the galaxies. We find that the background Lyman-limit flux due to galaxies might 
dominate (or be comparable to) the contribution from quasars at almost all redshifts 
if the absorption by HI in the ISM is neglected. The ISM HI absorption results in a 
severe diminishing of this flux — by almost three orders of magnitude at high redshifts 
to between one and two orders at z ~ 0. Though the resulting galaxy flux is completely 
negligible at high redshifts, it is comparable to the quasar flux at z ~ 0. 
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1 INTRODUCTION 

The study of the evolution of the ionization state of the dif- 
fuse intergalactic medium (IGM) is crucial to understanding 
the formation of structures in the universe. An important in- 
put in this study is the intensity and evolution of the back- 
ground hydrogen-ionizing flux. The universe is observed to 
be highly ionized at high redshifts (Gunn & Peterson f965, 
Giallongo et al. 1994) and it is widely believed that the 
main mechanism of ionization is photoionization by photons 
at and below the Lyman-limit. This view is borne out by 
estimates of background hydrogen-ionizing flux as inferred 
by the 'proximity effect' which fixes the value of this flux, 
J ion ~ i(r 21±0 - 5 ergs" 1 cm" 2 Hz" 1 sr" 1 (Bajtlik, Duncan & 
Ostriker f988; Bechtold f994; Giallongo et al. f996; Cooke 
et al. 1997) for 2 £ z Ss 4.5. It is not clear whether the 
observed quasars can provide sufficient contribution to the 
flux implied by this effect at z > 3 (Miralda-Escude & Os- 
triker f990; Giroux & Shapiro f996; Haardt & Madau 1996; 
Cooke et al. 1997). In light of this fact, it was speculated 
that high-redshift star-forming galaxies could be the dom- 
inant contributer to this background flux (Songaila, Cowie 
& Lilly f 990) . At smaller redshifts, there is a greater uncer- 



tainty in the level of this flux (Kulkarni & Fall f 993; Maloney 
f993; Vogel et al. f995; Reynolds et al. f995); however, a 



tight upper limit of Jion < 8 x f 23 erg s 



; Hz- 



at z ~ exists from Ha observations (Vogel et al. f995). If 
this flux is entirely due to the background hydrogen-ionizing 
photons then the contribution of quasars falls short of the 
observed value by nearly an order of magnitude at z ~ 
(Madau f992; Haardt & Madau 1996). 

Recent years have seen tremendous progress in under- 
standing the formation and evolution of galaxies at moder- 
ate and high redshifts. Prominent observations which made 
it possible are the probing of the evolution of galaxies with 
z < 1 by the Canada-France Redshift Survey (CFRS) (Lilly 
et al. f995), the discovery of high-redshift galaxies using 
the UV drop-out technique (Steidel et al. 1 996) in the Hub- 
ble Deep Field (HDF) (Madau et al. f996; Sawicki et al. 
f997; Connolly et al. f997), and a better understanding of 
the evolution of the neutral hydrogen content and average 
metallicity of DLA systems from z ~ to z ~ 4 (Lanzetta et 
al. f995, Pettini et al. f997). Several theoretical attempts 
have been made to explain these observations (Madau f 996, 
f997; Fall et al. f996; Frenk et al. f996; Guiderdoni et al. 
f997b, hereafter GHBM). These observations give valuable 
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clues about the history of star formation and the comoving 
cmissivity in several wavebands in the universe upto z ~ 4 
and have provided a unique opportunity to model and study 
the contribution of star-forming galaxies to the background 
hydrogen-ionizing flux (for recent attempts to estimate the 
hydrogen- ionizing flux from star-forming galaxies for z < 1.5 
see Giallongo et al. (1997); Deharveng et al. (1997)). 

The CFRS observations give the evolution of cosmic 
cmissivity at A ~ 2800 A (in the rest frame of the emitter) 
for z < 1 while the HDF observations estimate the emissiv- 
ity at A ~ 1620 A for 2.5 < z < 4. Our aim in this paper is 
to extrapolate the results from these important UV wave- 
lengths to calculate the emissivity (and background flux) 
at the Lyman-limit wavelength. The hydrogen-ionizing pho- 
tons are emitted by massive stars with lifetimes ~ 10 7 years, 
so their flux is not sensitive to the history of galaxy evolu- 
tion, which has much larger time scales. Conversely, one can 
view it as a direct measurement of the instantaneous SFR. 
The conversion of emissivity from the wavelengths probed 
by HDF and CFRS to 912 A requires knowledge of three 
factors: 

• the absorption of photons with A < 912 A within the 
stellar atmosphere. 

• the amount of dust and the spectral index of dust ab- 
sorption coefficient in the interstellar medium. 

• the amount of neutral hydrogen in the interstellar 
medium. 

The first factor is fixed by stellar atmosphere models and 
depends mainly on the choice of IMF which is guided by 
our ability to reproduce observed colours of galaxies. 

In extrapolating the amount of dust absorption from 
the observed value at 1620 A (or 2800 A) to the Lyman-limit 
wavelength, it is important to know the amount of dust in 
the galaxies as accurately as possible, because A912 1 Ay ~ 5 
as compared to Ai6oo/Ay — 2.5. Therefore a small uncer- 
tainty in the dust extinction at UV wavelengths translates 
into a much larger error in the extinction at the Lyman- 
limit. We can use FIR data available to fix the amount of 
dust at 2 ~ 0. At higher redshift, we can fix it by matching 
the observed average metallicity of DLA absorbers (Pettini 
et al. 1997) and by reproducing the cosmic infrared back- 
ground (CIB) reported by Puget et al. (1996). 

The most crucial factor in determining the absorption 
of hydrogen-ionizing photons in the ISM is the amount of HI 
and its distribution in galaxies. We can model this by using 
the available data on DLA absorbers which give the evolu- 
tion of the average HI content of these clouds for < z < 4 
(Lanzetta et al. 1995). In our study we have to assume that 
DLA absorbers belong to the same population of galaxies 
as the ones that are observed in emission in the CFRS and 
HDF. No clear consensus has been reached on this point yet, 
but recent observations may support the explanation that 
the DLA absorbers are either thick rotating disks or proto- 
galactic clumps (Haenelt et al. 1997; Prochaska & Wolfe 
1997). Thereafter We assume the HI, dust, and stars to 
be homogeneously distributed in the ISM. This assumption 
gives the maximum possible HI absorption and allows us 
to obtain a lower limit on the fraction of Lyman-limit pho- 
tons escaping from galaxies. We shall also briefly discuss the 
implications of the dumpiness of HI. 

We use a semi-analytic model of galaxy formation to 



study the evolution of background hydrogen-ionizing flux for 
~ 2 S3 4 (for details see, e.g. GHMB). In §2, we describe 
the relevant features this model and its various assumptions. 
The basic method of obtaining the background ionizing flux 
given the comoving emissivity and various sources of inter- 
galactic absorption is also briefly described in this section. 
Our results are detailed in §3. In §4 we discuss various un- 
certainties and possible improvements in our estimates and 
summarize our conclusions in §5. 



2 SEMI-ANALYTIC MODELS AND 
HYDROGEN-IONIZING FLUX 

We do not attempt to give a detailed description of semi- 
analytic models of galaxy formation and evolution. Readers 
interested in details are referred to the recent paper GHBM 
and references therein. We just introduce and define the pa- 
rameters required for our work. 

2.1 Collapse of halos and cooling 

We consider fluctuations in the matter density distribution 
at high redshift to be spherically symmetric ( "Top Hat" ap- 
proximation, see eg. Peebles 1980) on galactic scales. As a 
result of the violent relaxation process that they undergo af- 
ter collapsing at redshift z co ii, such halos virialize. One can 
then use the peaks formalism (Bardeen et al. 1986) to com- 
pute the mass distribution of such halos, as in Lacey and 
Silk (1991) and Lacey et al. (1993), for a given cosmologi- 
cal model. 

The baryonic gas subsequently radiates away its en- 
ergy, falling deeper and deeper in the potential well, until 
it reaches rotational equilibrium and settles into a disk-like 
structure. Thus, the length scale rdisk, of the cold gas ex- 
ponential disk can be inferred from the initial radius r v by 
conservation of angular momentum r d i S k — A,/r„ where Aj is 
the dimensionless spin parameter (Fall & Efstathiou 1980). 
Only disks can form in this formalism. The formation of el- 
liptical galaxies (and of bulges of spiral galaxies) has to be 
explained by the merging of disks. Kauffmann et al. (1993) 
and Cole et al. (1994) showed that this process can explain 
the current fraction of giant ellipticals among bright galax- 
ies. 



2.2 Star formation 

It is well known that local physical processes governing star 
formation are complex and depend on various parameters. 
Nevertheless, it seems that, on galaxy scales, one can phe- 
nomenologically define a mean Star Formation Rate (here- 
after SFR) per unit surface density, which is proportional 
to the total gas surface density (neutral plus molecular) di- 
vided by the dynamical time scale of the disk (Kennicutt 
1989, 1997). Hence we assume that the star formation time 
scale i* is proportional to the dynamical time scale of the 
disk: t-„ = (3(z)t dyn where t dyn = 2Tvr disk /V c . In the previous 
formula V c = (GM/r„)'' J is the circular velocity of the dark 
matter halo of mass M, and f3(z) is our first free parameter, 
that we allow to vary with redshift in order to reproduce 
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the shape of the cosmic comoving emissivity at 1620 A and 
2800 A from z = to z = 4. The SFR at time £ is then 
computed as follows: 



SFR(t) 



M gas (t) 
U 



(1) 



As shown in GHBM (their fig. 3), using such a prescrip- 
tion allows one to match the average value and width of the 
observed distribution of "Roberts time" for a z ~ bright 
disk galaxies sample given by Kennicutt et al. (1994), just 
by taking (3(0) ~ 100. Furthermore, as pointed out by these 
authors, the shape of the cosmic comoving star formation 
rate density curve (their fig. 9) suggests a strong evolution 
of f3 with redshift. Our model shares the same spirit: we 
consider the star formation history in the universe as a sum 
of contributions from a population of quiescent star-forming 
galaxies and a population of starburst galaxies whose respec- 
tive proportions are modelled through (3(z). We allow (3(z) 
to vary from ~ 0.5 to ~ 200 which respectively corresponds 
to starburst and quiescent mode of star formation. 



2.3 Stellar feedback 

Stellar feedback is modelled following the original work by 
Dekel and Silk (1986). It stems from massive star explo- 
sions which expel gas from the galaxies, preventing further 
star formation to occur. Physically, when the thermal energy 
ejected by supernovae becomes greater than the binding en- 
ergy of the halo, one expects such winds to be triggered with 
a certain efficiency. One can then compute the mass fraction 
of stars forming before the galactic wind: 

w (i+(v hot /v c rr\ (2) 



F* = 



(M» +M gas ) 



where Vhot and a are obtained through a fit based on SPH 
simulations of galaxy formation in which most of the feed- 
back effect is due to momentum exchange as in Cole et al. 
(1994). But one should bear in mind that there is much 
uncertainty on these parameters mostly due to the treat- 
ment of supernova remnant interactions with the interstellar 
medium. Moreover, as pointed out by Efstathiou (1992) and 
Blanchard et al. (1992), there is likely to be an overall high-z 
re-heating of the intergalactic medium which could prevent 
cooling in halos with circular velocities below V c ~ 20 to 
50 kms -1 , and possibly as high as ~ (200) 1 ^ 3 V C in case of 
adiabatic collapse. To sum up, the situation is very compli- 
cated because we lack a global theory of feedback processes. 
Therefore, we adopt a rather crude but simple approach to 
the problem: we do not attempt to model the redshift de- 
pendence of local or global feedback in this paper, but we 
just consider the values derived from numerical simulations 
for a typical ten percent feedback efficiency of momentum 
exchange. They give a = 5 and Vhot = 130 km s _1 . In the 
following, we keep the same value for a and allow our second 
free parameter Vhot to vary between 100 and 130 km s _1 . 



2.4 Stellar populations 

We use a coupled model of spectrophotometric and chemi- 
cal evolution in order to compute the age dependence and 
metallicity of the gas content, and of the UV to NIR spec- 
tra of the stellar populations in a self-consistent way. The 



end-product of these models are time dependent synthetic 
spectra and gas/metallicity evolution of galaxies. The model 
which is used here is based upon a new numerical scheme 
(isochrone), as well as on up-graded Geneva stellar tracks 
and yields, and will be described in Devriendt, Guiderdoni 
and Sadat (1998, in preparation). The photometric proper- 
ties of galaxies which are our main interest in this study are 
obtained after taking into account the intrinsic extinction 
due to both the neutral hydrogen absorption below 912 A 
and the dust (see following section). 



2.5 Dust and gas absorption 

Some of the energy released by stars in the UV and optical is 
absorbed by dust and HI gas (for wavelengths shorter than 
912 A ) and re-emitted in the IR and submm ranges. We 
would like to emphasize the difficulty of estimating these 
absorption processes since one has to address the compli- 
cated and crucial issue of the geometrical distribution of 
dust and gas relatively to stars. In the following, we assume 
that the gas is distributed in an exponential disk with mean 
HI column density at time t, 



<AWsm(£)> = M gas {t)/lAmHirr g 



(3) 



where r g = f c rduk is the truncation radius of the gaseous 
disk, f c being our third free parameter. The factor 1.4 ac- 
counts for the presence of helium. The mean optical thick- 
ness inside r g is given by: 



( Z g (t) \ ( (N m ,iSM(t)) 
v/z Q \ Z Q ) \2.\ 10 21 at cm- 
A ^ 3 



+ (AWsmW)^^-) 0(912 A), (4) 
where the first term contains: 

• the extinction curve for solar metallicity as measured 
by Mathis et al. (1983). 

• the gas metallicity (Z g (t)) dependence of the extinc- 
tion curve which is modelled as in Guiderdoni and Rocca- 
Volmerange (1987) and Franceschini et al. (1991, 1994), 
according to power-law interpolations based on the Solar 
Neighbourhood and the Magellanic Clouds, with s = 1.35 
for A < 2000 A and s = 1.6 for A > 2000 A. 

The second term is due to the hydrogen absorption and 
hence is only used for wavelengths shorter than 912 A (O 
being the Heavyside function), a H i = 6.3 x 10~ 18 cm 2 is the 
hydrogen ionization cross section at the threshold. 

Light scattering by dust grains is also modelled and 
the final optical depth is averaged over inclination angle as- 
suming stars, dust and gas are homogeneously mixed as in 
Dwek and Varosi (1996). These authors extended the ra- 
diative transfer calculation in a spherically symmetric ISM 
to the case of an ellipsoid which is more relevant for our 
study. The results are quite similar to those obtained with 
a classic "slab" geometry where stars, dust and gas are dis- 
tributed with the same scale height, and seem to satisfy the 
estimates of extinction given by Andreani and Franceschini 
(1996). As the galaxy is optically thick to Lyman- limit pho- 
tons, this homogeneous distribution of stars and dust implies 
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that only a tiny fraction (~ l/(rgi2)) of Lyman-limit pho- 
tons will be able to escape from the galaxy: the ones emitted 
by stars located in a thin outer shell (with column density 
~ 1/cthi — 1.6 x 10 17 at cm~ 2 ) of the gaseous disk. This 
procedure enables one to compute the bolometric luminos- 
ity absorbed by dust and re-processed in the IR/submm. 



2.6 Dust emission 

Here, we will not give details on how we derive the 
IR/submm spectrum (see GHBM), but just briefly outline 
the steps that lead to the emission spectrum. As mentioned 
in the previous section, we compute the optical depth of the 
disks, and the amount of bolometric luminosity absorbed 
by dust. The last step consists in the computation of emis- 
sion spectra of galaxies in the IR/submm wavelength bands. 
This is completed by using a three-component dust model 
(Polycyclic Aromatic Hydrocarbons, Very Small Grains and 
Big Grains), which is described in Desert et al. (1990). The 
method employed is very similar to that developed by Maffei 
(1994), and is based upon observational correlations of the 
IRAS flux ratios with the total IR luminosity (Devriendt, 
Guiderdoni & Sadat 1998, in preparation). 



2.7 Hydrogen-Ionizing Flux 

Once one has computed the fraction of hydrogen ioniz- 
ing photons that effectively escape from the galaxies, one 
still has to estimate the optical depth of the inter-galactic 
medium. This section describes the prescription we use. 
Given the comoving emissivity of hydrogen-ionizing flux 
e(u, z), the specific intensity of the background flux J(u , zo) 
(in ergs -1 cm~ 2 Hz -1 sr _1 ) is calculated using (Bechtold et 
al. 1987): 



T , , cH- 1 [°° (1 + 

J( ^ o) = IT/ dz (TTW 



x e(u, z)(l + z) 3 exp[-T(v ,zo,z)\ 



(5) 



Here t(vo,zo,z) is the optical depth for a photon emitted 
at redshift z with frequency v and observed at redshift zo 
with frequency vo = v(l + zo)/(l + z). In this paper we 
assume that this optical depth is due to HI in Lyman-a 
system clouds encountered along any line of sight. The ab- 
sorption due to singly ionized helium is negligible in deter- 
mining the hydrogen ionizing flux. Neutral helium is always 
a sub-dominant species in the Lyman-a clouds in case the 
background flux is dominated by quasars or star-forming 
galaxies. Therefore we neglect the absorption due to neutral 
helium. We also neglect the absorption due to diffuse HI in 
the intergalactic medium, because there is no evidence of its 
presence. 

The line-of-sight average optical depth r Ly (^o, zo, z) due 
to Poisson-distributed Lyman-a clouds can be expressed as 
(Paresce et al. 1980): 



TLy(fO, zq,z) 



pz poo 

/ / dzdNm,iGMV(Nm,iGM,z) 

J z JO 

x {1 - exp[-A?m,iGMO-HiW] }• (6) 
Here JVhi.igm is the neutral hydrogen column density of 



the Lyman-a clouds. V(Nhi,igm,z) is the average num- 
ber of clouds with HI column density between TVhi.igm and 
-/Vht.igm + d-/Vm,iGM in a redshift range z to z + dz along any 
line of sight. With the existing observations, it is not possible 
to uniquely fix the form of V(Nm,iGM, z) from z 4. 
For 2 > 1.7 we take model 3 of Giroux & Shapiro (1996) for 
Lyman-a clouds, with iV H i < 1.5 x 10 17 cm~ 2 . For z < 1.7, 
we take the redshift evolution given by Bahcall et al. (1993) 
and normalize the number of clouds per unit redshift at z = 
1.7 with model 3 of Giroux and Shapiro 1996. This is consis- 
tent with the findings of Bahcall et al. (1993) for clouds with 
equivalent width, W > 0.32 A {N m = 1.5 x 10 14 cm" 2 ). For 
Lyman-limit systems (TVhi > 1.5 x 10 17 cm~ 2 ), the results 
of Stengler-Larrea et al. (1995) are taken for < z < 4. 
The emissivity e(y, z) is determined from the galaxy evolu- 
tion models which fulfill a certain criterion described in the 
next section. The contribution of ionizing flux from 'repro- 
cessing' inside Lyman-a clouds is neglected in our estimates 
(for details see Haardt & Madau 1996). 



3 RESULTS 
3.1 Parameters 

In the following, we consider the so-called Standard Cold 
Dark Matter (SCDM) model with H = 50 km s _1 Mpc~\ 
fio — 1, A = 0, and as = 0.67. We take the baryonic fraction 
to be fig = 0.05, consistent with primordial nucleosynthe- 
sis. Our approach in this paper is to work within the frame- 
work of one cosmological model and a single set of values 
of the parameters, and to assume that all the changes that 
may occur from altering these choices can be compensated 
by appropriately adjusting the free parameters of the semi- 
analytic models of galaxy formation. As a matter of fact, it 
seems that the uncertainties from various galactic evolution 
processes like interactions, mergers, etc. are more important 
than the potential influence of the background cosmology 
and structure formation models (Heyl et al. 1995). 

For simplicity, we take a Salpeter Initial Mass Func- 
tion (hereafter IMF) with index x = 1.35. Stars have masses 
0.1 < m < 120M©, and we assume that the mass fraction 
forming dark objects with masses below 0.1 M a is negligible. 
Though we have used a Salpeter IMF throughout this paper, 
it is possible that the IMF at high redshift is top-heavy and 
hence produces more UV flux. We ran a few models with a 
top-heavy IMF at high-z and noticed that we overestimate 
the metal production by a considerable amount. However, 
these models may be allowed if a large fraction of metals are 
blown away with galactic winds in the IGM. Present obser- 
vations suggest that as much as half the metals at high red- 
shifts are in small column density Lyman-a clouds (Songaila 
& Cowie 1996), so it may not be unrealistic to assume that 
a substantial fraction of metals are blown from the galaxies. 
However, it can enhance the Lyman-limit flux by a factor 
of a few units at high redshifts which does not affect our 
results significantly (see the discussion the section 3.3). 

In order to put robust bounds on the hydrogen-ionizing 
emissivity t(y, z) coming from galaxies, we define two models 
which are representative of the uncertainties in the high-z 
measurement of the 1620 A cosmic emissivity (Madau et al. 
1996, Sawicki et al. 1997). Model I gives a fairly good fit 
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of Madau et aZ.'s results at high redshift. Model II matches 
values determined by Sawicki et al. (1997). Both models 
are completely denned by the three parameters previously 
mentioned. The first free parameter (the SFR efficiency pa- 
rameter), f3(z) is shown in Fig. |l| for both models. Its strong 
evolution with redshift is necessary to reproduce the sharp 
decrease of the cosmic comoving emissivity observed by the 
CFRS but the results are not very sensitive to the shape of 
the function used provided the peak is located at the right 
position 2 — 1.5 and that 0.5 < (3(z) < 200. For instance, 
Fig. |l| clearly shows that even though (3{z) evolves more 
strongly in model II (dashed curve) than in model I, there 
is no drastic difference in the cosmic emissivity around the 
peak in both models. 

The value of the concentration parameter f c is quite 
poorly known but, as pointed out in GHBM (their figure 6), 
the prescription given above allows one to match quite well 
the measured FIR and mm emission of a flux-limited sam- 
ple dominated by mild starbursts and luminous IR galax- 
ies (Franceschini and Andreani 1995) provided one takes 
f c ~ 2.7. However, such observations tend to pick out galax- 
ies with central, dense regions undergoing star formation 
(Sanders & Mirabel 1996) and therefore are biased to give 
smaller values of f c . On the other hand, as argued by Mo 
et al. (1997) and Lobo and Guiderdoni (1998) (in prepa- 
ration), semi-analytic models predict too small disk radii in 
the SCDM as compared to observations. Hence, we fix the 
value of f c to get the correct distribution of DLA (Fig. ^ 
which is based on the absorption properties of these clouds. 
/ c 's influence concerns mainly the HI distribution : the larger 
the value of this parameter, the more extended the HI ex- 
ponential disks and the fewer the DLA systems. There is 
likely to be a redshift dependence of this parameter, simply 
due to the fact that mergers have a well-known tendency to 
concentrate the gaseous content of disks. However, we did 
not try to model that effect but rather took an average value 
of the parameter over redshift. Averaged values f c = 7 for 
model I and f c = 5 for model II, allow us, as shown in Fig. 
to reproduce quite well the column density distribution of 
DLA at several redshifts. 

Finally, as discussed in the last section, the value of the 
third parameter (the feedback parameter), Vhot, is difficult 
to determine. However, it has very little influence on the 
results because it mainly affects the faint-end slope of the 
HI mass function, so we just tuned it in order to get the right 
amount of HI in DLA systems at high-z (Storrie-Lombardi 
et al. 1996). Such a requirement yields values for Model I 
and II of 120 and lOOkms -1 respectively. 

3.2 Luminosity Densities 

Figs. |§] and ^ show results for the two models discussed 
above. The "merit" of a model is judged by its ability to pre- 
dict correctly the evolution of emissivities in various wave- 
bands, the evolution of Qhi and the average metallicity of 
DLA absorbers from ~ z < 4. The predicted evolution 
of the Lyman-limit emissivities — with or without absorp- 
tion by HI in the interstellar medium of galaxies — is also 
shown in the figures. As seen in Figs. [] and ^, the two mod- 
els capture fairly well the broad features of the universe 
upto z ~ 4. The only discrepancy between the observations 
and theoretical estimates is that these models predict too 



high a value of Qui at z ~ (Figs. 3a and 4a). To inves- 
tigate this point in detail, we show in Fig. |^ the HI mass 
function at z = (Zwaan et al. 1997), along with the dis- 
tribution of HI column density in DLA absorbers at three 
redshifts (Lanzetta et al. 1995). Although the agreement at 
higher redshift is seen to be quite good, the disagreement at 
2 = can be noticed in Fig. 2d. Whereas the deviation of 
theoretical predictions from the observed HI mass function 
at small HI masses is probably due to our feedback pre- 
scription, we also overestimate slightly the HI content for 
1O 9 M < Mm < 5 x 1O 1O M . The disagreement at large 
HI masses is more important because much of the hydrogen- 
ionizing photons are seen to be emitted by the largest masses 
in our models. Nevertheless, this discrepancy, which would 
result in an overestimate of the HI absorption in the ISM in 
evaluating the hydrogen-ionizing emissivity, is somewhat al- 
leviated by the fact that the star formation rate is also higher 
because of its proportionality to the gas mass (Eq.([j])). 

The emissivity of Lyman-limit photons depends sensi- 
tively on dust and HI absorption in the ISM. At z ~ 0, the 
amount of dust in galaxies is fixed by matching the energy 
re-radiated in the FIR and mm bands with IRAS results at 
60 /im, as shown in Figs. |^ and|i|. At high redshifts, there is a 
greater uncertainty in the average amount of dust in galax- 
ies. We use Eq. (k|), which relates the dust optical depth 
with the average hydrogen column density (Nm,iSM(t)) and 
metallicity Z g (t) of a galaxy, to fix the dust content in a 
galaxy to within the observational uncertainties in the evo- 
lution of (./Vhi,ism(£)) and Z g (t). The recent detection of the 
CIB (Puget et al. 1996; Guiderdoni et al. 1997a) can also 
be used to constrain the history of dust emission in the uni- 
verse. GHBM used semi-analytic models of galaxy formation 
to explain the observed CIB. We plot results of their mod- 
els (A) and (E), which successfully explain Puget et al.'s 
results, in Figs. (3b) and (4b) along with the prediction of 
our models at A = 60 /im. It is quite clear that there is a 
good agreement between the two results. Therefore, we be- 
lieve that our method allows us to get a fair estimate of dust 
absorption for < z < 4. 

As is evident from Figs. |^ and ^, the HI absorption 
in the ISM is the most important factor in determining the 
hydrogen-ionizing emissivity. It is customary to compute the 
emissivity of galaxies and multiply by a constant escape frac- 
tion to account for this uncertainty (see e.g, Giallongo et al. 
1997). We model this by using Eq. rt4j) for the optical depth 
of the neutral hydrogen in a galaxy. As discussed above, in a 
geometry where the gas and stars are homogeneously mixed, 
it just means that a factor 1/(t912) of all the hydrogen- 
ionizing photons can escape the galaxy. As (T912) 2> 1, the 
validity of such an assumption depends crucially on the rel- 
ative distribution of stars and neutral hydrogen in the ISM. 
However, it should be pointed out that a homogeneous dis- 
tribution gives an overestimate of the absorption, i.e. an 
optically thick but clumpy medium will always result in less 
absorption, and therefore the average emissivity we estimate 
is a secure lower limit on the hydrogen-ionizing emissivity 
from star-forming galaxies. 

3.3 Galactic Background Lyman-limit Flux 

The evolution of the background flux of Lyman-limit pho- 
tons for various cases discussed above is shown in Fig. 
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along with the available observations at low and high red- 
shifts. The high redshift data shown in Fig. ^ is taken from 
Cooke et al. (1997). Other high redshift proximity effect 
calculations such as the one by Giallongo et al. (1996) give 
Ji on ~ (5 ± 1) x 10~ 22 erg s" 1 cm -2 Hz -1 sr" 1 independent 
of the redshift for 2 < z < 4. The evolution of the quasar 
contribution is also shown, based on the quasar luminos- 
ity function derived by Pei (1995). All the observational re- 
sults discussed before give the background hydrogen-ionizing 
flux, Jion, as opposed to the Lyman-limit flux we calculate. 
They are related as: Ji on = 3/(3 + 5) x J912, where <5 is the 
spectral index of the UV background flux (Miralda-Escude 
& Ostriker 1992). For quasar and galaxy UV backgrounds, 
S ~ 1-2 (see for instance Miralda-Escude & Ostriker 1990). 
We scale all the values of Ji on assuming S = 1, at all red- 
shifts, in Fig. ^ We also show the flux without dust absorp- 
tion to gauge the relative importance of dust absorption as 
compared to HI absorption. It should be pointed out that 
the curves without dust absorption are not normalized to 
HDF and CFRS fluxes at other wavebands so they cannot 
be taken as realistic estimates of Lyman-limit flux; their only 
purpose is to estimate the relative importance of dust and 
HI extinction of the Lyman-limit flux. As is clearly seen, 
dust absorption decreases the HI flux typically by a factor 
of a few while the HI absorption can diminish it by nearly 
three orders of magnitude at high redshift. 

If the absorption of Lyman-limit photons by HI in the 
ISM of a galaxy is neglected, then the background hydrogen- 
ionizing flux from star-forming galaxies can be comparable 
to or dominate the flux from quasars and might make up 
for the extra flux which might be required to explain the 
proximity effect (Cooke et al. 1997). However, we believe 
that it is highly unrealistic to completely neglect the HI ab- 
sorption in the ISM. As Fig. |s| shows, the ionizing flux at 
z ~ 3 decreases by more than three orders of magnitude 
when this absorption is taken into account. Such a decrease 
cannot be compensated by the uncertainties in our analysis 
and therefore it seems unlikely that star-forming galaxies 
dominate the hydrogen-ionizing flux at high redshifts. The 
situation is quite different at smaller redshifts. At z ~ 0, 
even the lower limit to galactic flux, (see the lower set of 
curves in Fig. |H|) might be more than the quasar contribu- 
tion. As discussed above, we overestimate the HI absorption 
in the ISM and so the value of the ionizing flux at z ~ 
could be higher by a factor of a few over the values shown in 
Fig. [| Therefore, though the quasars are likely to dominate 
the background hydrogen-ionizing flux at high-z, the sub- 
stantial contribution to this flux at smaller redshift is very 
likely to be due to star-forming galaxies. Our results are in 
qualitative agreement with the recent estimates of Giallongo 
et al. (1997) and Deharveng et al. (1997). 

The nature of the ionizing background at high red- 
shifts is also indicated by the recent observation of He II 
at high redshifts (see e.g., Davidsen et al. 1996). These ob- 
servations allow one to estimate of the softness parameter, 
Sl = J912 /J228 • At z ~ 3, Sl — 40 for quasar-dominated 
background flux (Haardt & Madau 1996) but can be as 
high as 1000 for background flux dominated by star-forming 
galaxies (Miralda-Escude & Ostriker 1990). However, var- 
ious uncertainties in the modelling of Lyman-a clouds at 
high- z as well as in the observations do not allow a firm 
conclusion on the softness of the ionizing background at 



high-z (Sethi & Nath 1997). We also note that the re- 
cent numerical simulations of Lyman-a absorbers show a 
good match with the observations if a low value of J912 
(~ 1-4 x 10 -22 erg s" 1 cm~ 2 Hz" 1 sr" 1 ) and a quasar-like 
spectrum is adopted (Miralda-Escude et al. 1996, Zhang et 
al. 1996). Therefore, it seems that apart from a few estimates 
of the proximity effect (Cooke et al. 1997), the value of ion- 
izing flux implied by the observed quasars may suffice to 
explain other observations at high-z, which is in agreement 
with our conclusions. 



4 DISCUSSION 

What are the various sources of uncertainties and errors in 
our estimates? We have assumed that the stars, dust, and 
HI are homogeneously distributed in the galaxy, which, as 
already discussed, is not entirely justified. In this section, wc 
discuss the possible implications of relaxing this assumption. 

Throughout this paper we use Eq. (|l| and (^) which give 
the SFR and optical depth averaged over the entire disk. We 
experimented with a few cases in which we considered the 
effect of local SFR and 7>. We subdivided the disk in shells 
of increasing radii and estimated the SFR and optical depth 
averaged over these shells; this enabled us to calculate the 
amount of Lyman-limit photons escaping from the disk as a 
function of the radial distance. In this case, both the SFR 
and the optical depth decrease towards the outer regions 
of the galaxy. The aim of this exercise was to investigate 
whether the decrease in SFR is less important than the de- 
crease in the optical depth, and a relatively larger fraction 
of Lyman-limit photons could escape from outer parts of a 
galaxy. However, this does not happen because, for our val- 
ues of f c , the disk remains optically thick to Lyman-limit 
photons even at the truncation radius r g where the HI col- 
umn density N aiilSM (r g ) = (-A/hi,ism)/1-6x/| exp(-1.6x/ c ), 
which means that even there the fraction of Lyman-limit 
photons which can escape remains ~ l/rgra- As both the 
SFR and T912 are proportional to (A'hi.ism) (Eqs. (|^) and 
(^)), the decrease in the optical depth is compensated by a 
similar decrease in SFR. Therefore, the total emergent flux 
from a galaxy remains the same whether we consider quanti- 
ties averaged over the entire disk or we sum the contributions 
of all the shells. 

Though it is extremely difficult to accurately model the 
dumpiness of the ISM, we give below qualitative argument 
to gauge its effect on the HI ISM absorption. Let us assume 
that all the HI in the ISM is in the form of HI clouds with 
number of clouds per column density /(JVhi.ism) distributed 
as: 

/(AWism) = , A f M = AN-;! ISU (7) 
al v hi, ism 

The optically thin clouds (JVhi.ism < ^h/) mimic the homo- 
geneous part of the ISM while the clouds with JVhi.ism > a hi 
correspond to the clumpy ISM. The total average optical 
depth, along any line of sight, due to Poisson-distributed 
clouds with this distribution is: 

/•• iv HI,ISM( max ) 
(t"ism) = / [1 - exp(— JVhI.ISMO'hi)] 

"'-?VHI,ISM(mm) 

X /(Ahi.ism) dAm.ISM. (8) 
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Our aim is to calculate the absorption from the clumpy 
part of the ISM and see what difference it would make if 
the same amount of matter was distributed homogeneously. 
To do that, we are interested in the ratio of the aver- 
age optical depth of the clumpy medium from Eq. (^) and 
(t^m) = {-/V H i,ism)o"hi, the optical depth if the HI in the 
clumpy medium was distributed homogeneously. We further 
assume that the ionizing stars are located outside the clouds 
in order for the HI absorption to be minimum. For the dis- 
tribution given by Eq. (M), we get: 



( r ISM ) 


(-7 + 2) 


'iVil+^Crnax) 


-N&£L(min)" 


1 


( r i^M ) 


(-7 + 1) 


"ivS+ 2 i(max) 


-N^+^mm) 


Cm 



(9) 



Observations of HI clouds in the ISM of our galaxy suggest 
the clouds have a distribution given by Eq. (bj) with 1.6 < 
7 < 2.2 for 10 18 at cm" 2 < TVhi.ism < 10 22 at cm -2 (Dickey 
& Garwood, 1990), which corresponds to the clumpy part 
of the ISM. For these values, 0.01 < (rf SM ) / (r" SM ) < 0.03 
(7 = 2 is excluded as Eq. ([]) is not valid for this value) . Our 
simple analysis suggests that the average optical depth from 
the clumpy part of the ISM is negligible as compared to the 
homogeneous part; and much of the absorption is caused by 
the fraction of HI in the homogeneous medium (or optically 
thin clouds) in the ISM. There is a great amount of un- 
certainty in determining the fraction of the homogeneously 
distributed HI in the ISM. If, for the purpose of this paper, 
we roughly assume that the HI is distributed equally be- 
tween the homogeneous and the clumpy ISM, as might be 
the case for our galaxy, then our analysis suggests that the 
average optical depth of the ISM is halved as compared to 
the case when the entire ISM is made up of homogeneously 
distributed HI. 

It is known that a significant fraction of young stars 
form inside optically thick clouds which suggests that the 
arguments given above underestimate the absorption. On 
the other hand, most of the ionizing stars form in OB asso- 
ciations. The UV photons from these stars could puncture 
the layer of neutral hydrogen to escape into the galactic halo. 
This fact can be used to explain the existence of HII at large 
scales heights (Reynolds 1984) which requires a large frac- 
tion of the photons to escape in the halo (Dove & Shull 
1994). However, it is not clear what fraction of this ioniz- 
ing flux can escape the halo of a galaxy. The observational 
status of hydrogen ionizing flux from other galaxies is uncer- 
tain. Recent observation of 4 low redshift starburst galaxies 
(Leitherer et al. 1995) allows one to get only upper limits 
of 57%, 5.2%, 11.3%, and 3.2% on the escape fraction (Hur- 
witz et al. 1997). Other arguments based on Ha luminosity 
density of the universe give a more stringent upper limit of 
~ 1% on the escape fraction (Deharveng et al. 1997). We 
note that our simple assumption, without puncturing, al- 
ready gives an average escape fraction of 0.4-1% at small 
redshifts, which is within the observational uncertainties. 
Therefore, this effect seems unlikely to be more than a few 
percents at low redshift. However, at high-2, as the star for- 
mation rate is higher, it might well be that the puncturing 
effect allows a larger fraction of ionizing photons to escape. 

Throughout this paper we assumed that the damped 
Lyman-a systems are progenitors of present day galaxies. 
This assumption has recently been questioned by Prochaska 



and Wolfe (1997). They studied the kinematical properties 
of damped Lyman-a systems at high redshift and concluded 
that most damped Lyman-a systems correspond to thick 
rotating disks with rotational velocities ~ 225 km s . We 
compare these results with the predications of our model. 
We notice that for standard CDM model, even though about 
ten percent of the objects have rotational velocities beyond 
the value of ~ 225kms _1 , the average rotational velocity 
of damped Lyman-a systems is 120kms _1 at 2 ~ 3, which 
is too small to be consistent with Prochaska and Wolfe's 
results. However, this discrepancy is not serious for vari- 
ous reasons. First, the relationship between rotational ve- 
locity and iVm depends sensitively on the cosmology. We 
tried a structure formation model with an open universe 
which has £1 = 0.3. In this case, the average rotational ve- 
locity increases to 160 kins -1 with more than one fourth of 
the objects with velocities greater than ~ 225 km s -1 . This 
increase is not unexpected because in open universes struc- 
tures of galaxy sizes collapse at higher redshifts as compared 
to flat cosmology. Second, the distribution of rotational ve- 
locities for damped Lyman-a systems depends crucially on 
the feedback. An increase in feedback parameter Vhot results 
in an increase in the average rotational velocity by cutting 
off contributions from smaller systems. Such an increase in 
the feedback might come from either local uncertainties in 
feedback (Eq. (2)) or from a global process like a hot IGM 
with temperature T ~ 10 6 K at high redshifts and therefore 
would depend on the thermal history of the IGM. We also 
compare the distribution of impact parameters in our mod- 
els with the inference of Prochaska and Wolfe and found 
good agreement for both cosmogonies. Also, other recent 
analyses suggest that the evidence of damped Lyman-a sys- 
tems being thick rotating disk is not compelling (Ledoux 
et al. 1997a, Ledoux et al. 1997b.). Ledoux et al. stud- 
ied the kinematical properties of a sample with 26 damped 
Lyman-a systems and concluded that whereas the velocities 
upto 120kmsec _1 might correspond to rotations of individ- 
ual systems, higher velocities probably involve more than 
one component (Ledoux et al. 1997a, Ledoux et al. 1997b). 
A similar conclusion, based on simulations, was reached by 
Haehnelt et al. 1997. Therefore, it seems that though the 
recent studies have thrown light on the nature of damped 
Lyman-a systems at high-2, it is still too early to draw 
firm conclusions. However, should it turn out that damped 
Lyman-a systems at high 2 correspond to merging proto- 
galactic HI clumps rather than rotating disk, our conclusions 
about the HI absorption inside the ISM of a high-2 galaxy 
and consequently the Hydrogen-ionizing flux will be signifi- 
cantly weakened in the sense that it will be more difficult to 
support the view that sites where a high HI column density 
is detected also correspond to star forming regions. 



5 CONCLUSION 

We estimated the contribution of star-forming galaxies to 
the background Lyman-limit flux taking into account the 
HI and dust absorption in the ISM of individual galaxies 
in a self-consistent way with the cosmic star-formation his- 
tory. We assumed that DLA systems correspond to star- 
forming regions at high-redshift. We conclude that while 
star-forming galaxies are unlikely to dominate the back- 
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ground hydrogen-ionizing flux at high redshift, they are most 
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likely 1 ) do so in the present universe. The current uncer- 
tainties of modelling do not allow us to calculate the redshift 
of cross-over from quasar-dominated to galaxy-dominated 
background flux. As already discussed, a good discrimina- 
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tor bet 



tweeu Lhtise Lwu suuices uf iunizhig-bdLkgiuuiid is the 
softnes a of their apectra. Future high rcaolution obacrvationa 

of metals (e.g. Carbon) in their various ionization states in 
the low-column density Lyman-a clouds for ~ z < 4 might 
enable one to find the transition from quasar-dominated to 
galaxy-dominated background flux. 
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Figure 1. Efficiency parameter evolution with redshift. The thick solid and dashed lines show the values used for model I and II 
respectively. Analytical expressions for (3(z) are 0(z) = 50(1 + 1(T 6 (1 + z coll ) 13 )/(l + z coll f (model I) and f3(z) = 10(1 + 10" 5 (1 + 
z C ou) ie ')/{l + z coll f (model II). 
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Figure 2. Fig. 2a: Distribution of damped Lyman-o absorbers as a function of their HI column density at z ~ 1. Thick solid and 
long-dashed lines are respectively the predictions of models I and II. Fig. 2b: Same as panel a) for z ~ 2 . Fig. 2c: Same as panel a) for 
z ~ 3. Fig. 2d: HI mass function of galaxies at z ~ 0. Thick solid , long dashed , and dotted curves are respectively models I and II 
predictions and a Schechter fit to the data given by Zwaan et al. 97. 
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Figure 3. Fig. 3a: Evolution of the cold gas density parameter in DLA absorbers. The fit from Model I is indicated by the thick solid 
line. Fig. 3b: Rest— frame comoving luminosity density from optical to FIR. The thin dotted , dot-dashed , triple- dot- dashed and solid 
lines respectively stand for emissivities at 4400 A, 10000 A 22000 A and 60 (im given by model I. Open squares: local and Canada- 
France Rcdshift Survey (Lilly et al. 96). Open circle: 60 (im local density corresponding to one third of the bolomctric light radiated 
in the IR (Saunders et al. 90). Solid circle : 2.2 fim local luminosity density (Gardner et al. 96). The thick solid and dashed lines are 
respectively models (A) and (E) from GHBM that give reasonable fits of the CIB. Fig. 3c : Rest-frame comoving luminosity density 
from far— UV to UV. The thin solid and long-dashed lines, represent emissivities at 1620 A and 2800 A given by model I. The thick 
solid and long-dashed lines are predictions of the same model for luminosity at 912 A respectively without and with HI absorption, 
but with dust absorption. Solid diamonds: photometric redshifts in the Hubble Deep Field (Connolly ct al. 97) taking into account NIR 
data. Solid squares: other estimates of photometric redshifts in the HDF (Sawicki et al. 97). Open triangles: HDF with redshifts from 
Lyman-continuum drop-outs (Madau ct al. 96). Fig. 3d : Evolution of the mean metallicity in damped Lyman-a absorbers. Prediction 
from Model I is indicated by the thick solid line. Since the chemical evolution model is a closed-box one, the metallicity of the systems 
is certainly overestimated if a fraction of the metals is ejected in the IGM. 
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Figure 4. Same as Fig. 3 for model II. 
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Figure 5. The evolution of Lyman-limit flux is shown for the various cases discussed in the text. The solid lines give the predictions 
of Model I. The lines of increasing thickness correspond respectively to models without dust and HI absorption, with dust but without 
HI absorption, and with both dust and HI absorption. The dashed lines are the corresponding predictions of Model II. The dotted line 
shows the evolution of background Lyman-limit flux from quasars. The upper limit at z ~ is from Vogcl ct al. (1995), the data point 
at z ~ 0.5 is from Kulkarni and Fall (1993), the high redshift data points are from Cooke et al. (1997). 
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